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The Analysis of Global Cloud and Radiation Data for the 
Study of Cloud-Climate Interactions 

NASA Grant NAGW - 3150 

During the past year, graduate student Kathryn Ginger has analyzed ISCCP Cl data for 
two regions, the stratocumulus area off California and the tropical Atlantic. Her work has led to 
an M.S. thesis which is in the process of being submitted in its final form. Copies will be mailed 
to the Technical Officer and the NASA Center for Aerospace Information when it is ready for 
dissemination. 

The first part of the study dealing with stratocumulus clouds is complete and the results 
form the technical portion of this report. A manuscript that has recendy been prepared covering 
this subject is attached as an Appendix. The abstract describes the key finding that cloudy mean 
LWP is invariant with cloud fraction for cloud cover between 20% - 80%. With reference to die 
Appendix, this is shown in Figures 7a and 7b which are based on LANDS AT data. Similar 
conclusions are reached with ISCCP data shown in Figures 9 and 10. This has important 
implications for General Circulation Models. A possible application is described in the Discussion 
section of the Appendix. We can basically use an idea that has been proposed for dealing with 
GCM hydrology. Efforts are being initiated along this direction with the Goddard GCM 
developed by Dr. Yogesh Sud. We have collaborated with Dr. Bruce Wielicki of NASA Langley 
in this effort. The material will be presented at the Eighth AMS Radiation Conference in 
Nashville, TN in January 1994 and is expected to be published during calendar year 1994. 

During the next grant year, results from the latter part of Ms. Ginger's thesis will be 
prepared for publication. A preview of the results may be seen in the thesis which is being mailed 
separately. As per our original research proposal, a model of the small scale liquid water path 
distribution will be solved for the radiation field which will be compared with ERBE data which is 
at coarser resolution. We also intend to corroborate the LWP distributions found in our study 
with SSM/I data which will be made available to us by Dr. Grant Petty of our department. The 
level of effort that was requested originally has not changed. 

Graduate students supported: 

Kathryn Ginger (M.S., Dec. 1993) 

Jaya R. Rao (joined group in June 1993) 
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Abstract 


A study has been made of the reladonship between mean cloud radiative properties and 
cloud fraction in stratocumulus cloud systems. The analysis is of several LANDSAT images and 
three-hourly ISCCP C- 1 data during daylight hours for two grid boxes covering an area typical of 
a general circulation model (GCM) grid increment. Cloud properties were inferred from the 
LANDSAT images using two thresholds and several pixel resolutions ranging from roughly 1/16 
km to 8 km. At the finest resolution, the analysis shows that mean cloud optical depth (or liquid 
water path) increases somewhat with increasing cloud fraction upto 20% cloud coverage. More 
striking, however, is the lack of correlation between the two quantities for cloud fractions 
between roughly 0.2 and 0.9. When the scene is essentially overcast, the mean cloud optical 
depth tends to be higher. Coarse resolution LANDSAT analysis and the ISCCP 8 km data show 
lack of correlation between mean cloud optical depth and cloud fraction for coverage less than 
about 95%. 

This study shows that there is perhaps a local mean liquid water path (LWP) associated 
with partly cloudy areas of O (100 km). A method has been suggested to use this property to 
obtain the cloud fraction in a GCM when the model computes a grid box mean LWP. 
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Introduction 


The role of clouds in modifying the earth's radiation budget has never been doubted. In 
recent years, there has been the realization that climate simulations will not gain widespread 
acceptability as long as there are nagging doubts related to the manner in which clouds are treated 
in these models. An appreciation of the current status may be gained from the comprehensive 
study made by Cess et al. (1990) who attributed most of the differences in the radiative response 
of nineteen general circulation models to differences in their cloud parameterization. These 
models had been forced with an imposed positive and negative uniform change in sea surface 
temperature, a forcing that acted as a surrogate for the anticipated radiative forcing produced by 
changes in the concentration of radiatively active atmospheric constituents. 

Concern with this weak element in present climate models has led to the identification of 
the role of clouds in atmospheric and hydrologic systems as the area of scientific study of the 
highest priority (Committee on Earth Sciences, 1989). Current research includes the cataloging 
and interpretation of global cloudiness data in the International Satellite Cloud Climatology 
Project (ISCCP, Rossow and Schiffer, 1991) and parallel field observations (Cox et al., 1987). It 
is felt that the analysis of data from global and regional studies will help to formulate improved 
cloud generation and radiative parameterization schemes in numerical atmospheric models. 
However, in order to assimilate observational results, it is first necessary to set up a framework 
for the interpretation of cloud data presented to modelers. This is necessary because 
observational results are obtained through an inversion process based on several assumptions, 
some quite reasonable, but others primarily for the sake of expedience. Models, on the other 
hand, compute radiative fields directly, using given constituent properties and thermodynamic 
fields under the assumption that the atmosphere is stratified into horizontally homogeneous layers 
occupying a defined area. 

This study is devoted to the interpretation of information related to cloud fields that can be 
directly used in numerical models. We have restricted ourselves to selected regions and cloud 
formations and to only two sources of information. These are the analysis of LANDSAT images 


of stratocumulus (Wielicki and Welch, 1986) and ISCCP data from the stratus area off the west 
coast of North America. The area of the LANDS AT images analyzed is 58 km on a side, 
comparable to the grid increment of three dimensional mesoscale models while ISCCP cloud 
statistics are for a 2.5° x 2.5° latitude-longitude area which is typical of the spatial scale of global 
general circulation models (Giorgi and Meams, 1991). Further, we are confining ourselves to day 
time observations of visible reflectance since these provide a first order estimate of optical 
thickness. Even this limited sample reveals the necessity of establishing some ground rules for the 
meaningful application of satellite cloud data to atmospheric radiation models. Section 2 
introduces the nomenclature used in the study, in particular a definition of cloud fraction that is 
appropriate for the direct computation of radiation fields in climate models. Section 3 is an 
analysis of an ensemble of LANDSAT images of stratocumulus clouds within the framework 
defined in Section 2. Section 4 extends the analysis using ISCCP data for stratocumulus cloud 
fields. Section 5 summarizes the results and discusses the implications for cloud modeling. 

2. Cloud Properties 

Inspection of visible images of the earth from space show that cloudiness is ubiquitous. 
Cloudy areas stand out in contrast to the generally darker background, particularly over oceanic 
areas. Images constructed from thermal emission in the atmospheric window also reveal 
cloudiness through the temperature difference between cloud tops and the usually warmer surface. 
Quantitative analysis generally begins with a procedure to discriminate cloudy areas from the 
background clear areas. The fraction of an arbitrary horizontal area determined to be cloud 
covered is the cloud fraction. Further processing can provide information on the optical 
properties of the clouds identified in the prior step. Unfortunately there are no unambiguous and 
universal rules that can be followed to identify cloudy areas and ascribe optical properties to them. 
This is a consequence of the rich horizontal variability of cloudiness at all scales compounded by 
the presence of geometric structure. 

Nevertheless, in order to transfer information obtained from the analysis of remotely 
sensed observations into a cloud-radiation scheme to be used in an atmospheric model, it is 
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necessary to create a discrete model of cloudiness. This may be appreciated by considering the 
manner in which the presence of clouds is determined in current three dimensional atmospheric 
models. Generally, two classes of clouds are formed: stratiform clouds and convective clouds. 
The former are generated by a super saturation criterion in the grid box and the latter through a 
convective parameterization. The relative humidity threshold for the presence of super saturated 
clouds could be 100% (Randall et al., 1989) or some lesser value which is incorporated in a 
fractional cloud cover algorithm (Slingo, 1987). There are numerous different schemes for the 
diagnosis of convective clouds. In addition to the determination of the presence of clouds, 
models are now assigning optical properties to the generated clouds usually in a diagnostic sense 
using empirical relations involving temperature, cloud water or altitude. A tabulation of cloud 
generation techniques and optical thickness schemes incorporated in current general circulation 
models (GCMs) is provided in Cess et al. (1990). 

The radiative fluxes within the atmosphere and net flux at the boundaries for each model 
grid point are computed from the cloud cover, cloud optical properties and the temperature and 
gaseous constituent profile. The information necessary to compute radiadve fluxes is provided 
only at a spatial scale corresponding to the grid increment of the model, Ax, which is typically 
100-500 km. The atmospheric column within these boundaries is allowed to be clear or covered 
by extensive clouds over a prescribed fractional area (including complete overcast corresponding 
to a cloud fraction of 1.0). In the model envisaged here, the horizontal area representing the 
numerical model grid is considered to be subdivided into areas of side, Ap, which will be called 
the pixel size or resolution. It should be noted in passing that the grid increment of GCMs is also 
referred to as the model resolution, perhaps erroneously (Pielke, 1991). The distinction between 
pixel resolution and model grid increment should be kept in mind as it is central to the framework 
of this study. 

We now state that a pixel is allowed to be clear or cloudy but not partly cloudy, and if 
cloudy, consists of a column of non-zero optical thickness. The model grid of side Ax then 
contains several pixels which are clear or cloudy. The ratio of cloudy to total number of pixels is 
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defined as the cloud fraction of the model grid. It will be further assumed that the individual 
pixels contribute independently to the radiation flux at any vertical level of the model grid. 
Cahalan et al. (1993) have shown that this is a justifiable assumption for marine boundary layer 
clouds. The definition of cloud fraction presented here is energy conserving at any pixel 
resolution and more importantly direcdy applicable to GCM radiation computations. 
Unfortunately, attempts at inverting the reflected radiance field to fit this scheme are confounded 
by the nature of cloudiness (Schertzer and Lovejoy 1987; Cahalan and Joseph, 1989) which 
exhibits variability at all length scales. As a consequence, the pixel size and reflectance threshold 
above background used to discriminate clouds will determine the cloud fraction and pixel 
reflectances of a particular image (Wielicki and Welch, 1986; Wielicki and Parker, 1992; Chang 
and Coakley, 1993). However, the distribution of reflectances when averaged over the image and 
the contribution of the areas designated clear provides the grid mean radiative property required 
by modelers. 

A further consideration is the difference in temporal scale between satellite images, which 
provide essentially instantaneous fields of cloudiness and GCMs in which cloud fields are updated 
at intervals ranging from one to twelve hours. In some models, the radiation field is held fixed 
within these time intervals even though thermodynamic fields are allowed to change. The time 
scale is much longer than that associated with turbulent processes within clouds and even the life 
cycle of individual cells embedded within cloud fields. Again, what is required is a statistical 
description of the radiative properties of the pixels contained within the GCM grid averaged over 
the time interval between cloud diagnosis updates. 

Once we have established the ground rules stated above, a statistical description of 
cloudiness can be formulated. The usual statistical model of cloudiness describes the frequency of 
occurrence of various classes of sky cover and is built up from conventional meteorological 
observations. These models can also be extended to apply to horizontal scales that are larger or 
smaller than the observations (Falls, 1974; Henderson-Sellers and McGuffie, 1991). Application 
of these cloud models to global sampling strategies has focused primarily on cloud amount 


4 



(Hughes and Henderson-Sellers, 1983; Rossow, 1989), which by itself is not sufficient to describe 
the radiative properties of the area under consideration. Now that global estimates of cloud 
optical properties are available from ISCCP (Rossow and Schiffer, 1991), a model of cloudiness 
more appropriate for GCM applications is possible. 

A particularly simple model analogous to those used in hydrology (Entekhabi and 
Eagleson, 1989) considers the GCM grid area average radiative property, say the planar albedo, 
A, to be the expectation of that property, E(A). The model is illustrated for the total atmospheric 
column but in principle can be applied to a particular layer as well. Following Stephens (1988) 
and Rossow (1989), it will be assumed that the atmospheric column contained in the horizontal 
area (Ax) 2 as viewed from the top of the atmosphere is composed of variable column amounts of 
liquid water and also a region where the presence of condensate can not be detected by the 
applied threshold criterion (the clear region). For simplicity, let us assume that the clear area has 
zero albedo and that the cloudy portion consists of several pixels of area (Ap) 2 whose individual 
radiative properties can be measured from a remote sensing system. If Q Watts of radiant energy 
is reflected by the entire grid of area (Ax) 2 , then 


E(A) = 


Q 

M-o^o (Ax) 


( 1 ) 


where |i 0 is the solar zenith angle, S 0 is the solar constant and E(A) is the expectation or grid area 


averaged planar albedo. If a fraction, C, of the total number of pixels is cloudy, then the mean 
albedo of the cloudy pixels is E(A)/C and the pixel albedos may be represented by a distribution 
f(A). For example, the study of fair weather cumulus clouds by Wielicki and Welch (1986) 
suggests an exponential distribution of the form 


f(A)= 


C 


E(A) 


exp 


C A 
E(A)_ ’ 


A>0. 


( 2 ) 


When (2) is to be applied to a satellite image it is understood that A should be greater than some 
non-zero threshold for discriminating clouds from the background and that C will depend on this 
threshold. In general, for an arbitrary area (Ax) 2 , the quantities C, E(A) and the distribution f(A) 
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will vary from one instant to another and at a particular time these quantities will vary from 
location to location depending on the cloud field. Moreover, these quantities will also be different 
when the averaging area (Ax) 2 is changed or the pixel size, Ap, is changed. In this study the 
radiative properties considered are the nadir reflectance measured by the LANDSAT family of 
satellites, and the effective spherical albedo from ISCCP. 

In the analyses that will follow, cloud fraction and area mean cloud radiative (and optical) 
properties will be obtained using different thresholds and pixel resolutions. These will be 
compared to 'reference' properties which are obtained using the bispectral threshold and the 
highest resolution available. The differences are at times referred to as 'errors' but it should be 
understood that the 'reference' measurements are not absolute but are defined within the 
constraints of the assumption that the pixel is completely full of cloud of constant optical depth. 
Even for the highest resolution data available this is not strictly true. 

3. LANDSAT Analysis 

Wielicki and Parker (1992), hereafter denoted WP, showed that threshold cloud retrieval 
methods similar to ISCCP are subject to two primary error sources: an underestimation of the 
amount of optically thin cloud, which cannot be detected by the threshold, and an overestimation 
of cloud amount caused by triggering the threshold with partially cloud filled fields of view. The 
former effect dominates for cirrus and both are important for oceanic boundary layer cloud (WP). 

The study by WP, however, did not consider derived cloud reflectance or cloud optical 
depth. The consequences of using different thresholds and resolutions are summarized 
schematically in Figure 1 which is illustrated for measurements of nadir reflectance. It may be 
noted that the two lowest panels show the same cloud fraction ( 25 %) but the cloud reflectance 
and inferred optical depth at high resolution may exceed the low resolution value by a 
considerable amount. These differences will tend to increase as the regional cloud fraction 
decreases. Because of these uncertainties, we will use the high spatial resolution LANDSAT data 
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to check the accuracy of simulated ISCCP results for cloud reflectance and optical depth and to 
corroborate any conclusions drawn from ISCCP results. 

The LANDSAT Thematic Mapper (TM) instrument flew on the LANDSAT 4 and 
LANDSAT 5 spacecraft. The TM has seven narrow spectral bands with central wavelengths of 
0.48, 0.56, 0.66, 0.83, 1.65, 2.21, and 11.5 Jim. The six solar reflectance bands have a field of 
view of 28.5 meters, while the 1 1.5-jim band has a field of view of 1 14 meters. All seven bands 
are atmospheric windows (i.e., small absorption by atmospheric gases). 

As in WP, we selected the 0.83-|im and 11.5-tim spectral bands as similar to those 
available on meteorological satellite instruments such as the VAS (VISSR Atmospheric Sounder) 
imager on the GOES (Geostationary Operational Environmental Satellite) satellites and the 
AVHRR (Advanced Very High Resolution Radiometer) imager on the polar orbidng satellites. 
These near-visible and infrared window channels are similar to those used for the majority of 
satellite cloud retrieval algorithms, including the ISCCP (Coakley and Bretheron, 1982; Minnis et 
al., 1987; Rossow et al., 1985). While the 0.66-(im band would be a closer match to the GOES 
visible spectral band, the 0.83-pm band was selected because of its greater dynamic range and 
lower ocean surface reflectance. 

The reference cloud cover defined for this study is determined using a bispectral threshold 
method similar to ISCCP on 57 meter resolution LANDSAT data. For convenience, 57-meter 
resolution data will henceforth be referred to as 1/16 km resolution data. The 114-meter 
resolution 11.5-jim band data are replicated to provide bispectral data at 57-meter spatial 
resolution. The replication of the 11.5-fim data has little or no impact on the results since the 
clouds most sensitive to spatial resolution are found to be cumulus and stratocumulus. These 
boundary layer clouds are detected almost exclusively by the 0.83-p.m channel. The reference 
cloud cover uses a nadir reflectance threshold R, = R clr + AR and a brightness temperature 
threshold T, = T dr - AT, where R clr is the nadir clear-sky reflectance and T clr is the clear-sky 
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brightness temperature. Nadir reflectance is calculated as an equivalent Lambertian reflectance 
(Wielicki and Welch, 1986), 


where I is the nadir radiance. With this definition, a perfect Lambert reflector would have a nadir 
reflectance of 1.0. Calibration coefficients for the LANDSAT spectral bands are taken from 
Salomonson and Barker (1987) and Markham and Barker (1986). The actual reference 
thresholds vary as the background variability of the clear ocean increases or decreases. For the 45 
cloud fields used in the present study, AR has an average value of 0.017 and a standard deviation 
of 0.008. Corresponding values for AT are -1.5K and 0.6K. These values are within 0.001 and 
0.1K of those found in WP. Estimated rms error in the reference cloud fraction is less than 0.05, 
and an error analysis can be found in WP. 

The location and time of observation of the 45 scenes is listed in Table 1. The selected 
scenes are typical of the variety of boundary layer clouds expected in the subtropics, ranging from 
trade cumulus to solid stratocumulus decks. They were selected from browse images to obtain a 
wide range of cellular structures: open/closed, cloud cell diameters from 0.5 to 50 km, and cloud 
streets. After applying the LANDSAT threshold (LS) and the somewhat more stringent ISCCP 
threshold (IS) separately on each scene, the cloud fraction and cloud pixel mean nadir reflectance 
have been computed and are shown as a scatter plot in Figs. 2(a) and (b). In addition, the data in 
Fig. 2(a) are listed in Table 1 so that each point on the plot may be identified with the particular 
scene. As in the ISCCP Cl data, the average is over the cloudy pixels alone such that, with 

reference to Figure 1, the mean is of I,, I 2 , for each pixel identified as cloudy by the particular 

threshold. The general pattern is not influenced greatly by choice of threshold. For cloud 
fractions less than 0.2, the mean nadir reflectance increases with cloud fraction. There is then a 
wide range of cloud fraction over which there appears to be little correlation between reflectance 
and cloud fraction. Completely cloud filled scenes, however, tend to have the highest values of 
mean reflectance although there are exceptions to this rule. The most obvious is scene 11 which 
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is roughly 40% cloudy but has highly reflective clouds. This is a scene of cloud cover following a 
cold air outbreak in which the larger LANDSAT image shows essentially clear skies over half the 
image and a very thick unbroken deck over the rest of the scene. 

When the same scenes are analyzed at progressively coarser resolution but using the same 
threshold, there is a marked evolution in the pattern. This is shown in Figs 3(a) - (d) which shows 
a series of scatter plots of analyses at 1/16, 1/2, 2 and 8-km resolution respectively using the 
ISCCP (IS) threshold. Figure 2(b) and Figure 3(a) are identical with the other panels showing the 
effect of coarsening the resolution. Figure 3(d), which may be considered to be simulated ISCCP 
data shows that mean nadir reflectance and cloud fraction are uncorrelated at least out to a cloud 
fraction of 0.8. 

The effect of coarsening the resolution and changing from LS to IS threshold is 
summarized in Figures 4(a) and (b). As mentioned in WP, the coarser resolution analysis 
overestimates cloud fraction except for scenes with small amounts of scattered clouds. For these 
cases, the resolution results in the larger pixels being classified as clear. An extreme example is 
scene 17 which has a cloud fraction of 0.1 1 at 1/16-km resolution and using the LS threshold but 
is classified as completely clear at the 8-km resolution using the ISCCP threshold. However, it is 
important to realize that changes in cloud fraction are accompanied by concomitant changes in 
mean nadir reflectance. 

In order to examine the implications of the patterns exhibited in the panels of Figs. 2 and 3 
it is necessary to present the data in terms of a physical quantity such as the optical depth or cloud 
liquid water path. This is accomplished by the following procedure. For each cloud pixel, the 
small contribution of the ocean surface reflectance is first removed following the method of Platt 
et al. (1980). The remaining cloud reflectance is then compared to a look-up table to convert 
nadir reflectance into an estimate of cloud optical depth at 0.83 (im. The look up table is based 
on radiative transfer calculations using a multiple scattering finite difference model (Barkstrom, 
1976; Suttles, 1981), which has been checked for consistency against the adding doubling 
approach used by ISCCP. Single scattering properties are derived using Mie calculations tor a 
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size distribution of water droplets with an effective radius r e of 10 p.m and a variance of 0.10, 
consistent with ISCCP (Rossow et al., 1991). Figure 5 shows the resulting reflectance at nadir as 
a function of cloud optical depth for solar zenith angles of 25, 45 and 60 degrees, covering the 
range of values used in this study. A more complete description of the radiative calculations can 
be found in Wielicki et al. (1990). Finally, an estimate of cloud liquid water path (LWP) is made 
by using the following relationship (Stephens, 1984), 

LWP = - tr (4) 

3 e ' 

where x is the optical depth and r e is the effective particle radius. 

A linear average of the pixel optical depth or LWP of the cloudy pixels alone is then a 
measure of the amount of liquid in the region. This process was First carried out at the highest 
resolution using the LS threshold. After estimating the above parameters for the LANDSAT full 
resolution data we then spatially averaged the LANDSAT radiance data to a spatial resolution of 
7.2 km, approximating the ISCCP GOES data source. As in WP, we then applied the ISCCP 
cloud detection thresholds to the spatially averaged data, and computed ISCCP estimates of cloud 
reflectance, cloud optical depth, and cloud LWP. 

As mentioned earlier, estimates of mean reflectance are influenced by resolution and 
threshold. This of course carries over to estimates of the optical depth. Figure 6 is a comparison 
of the linearly averaged optical depth for the 45 scenes at two resolution and using the LS and IS 
thresholds. The coarse resolution analysis underestimates the optical depth for partly cloud 
scenes. This underestimate could be substantial. 

The inferred mean LWP for the scenes is presented in Fig. 7(a) and (b) as a function of the 
regional cloud fraction. The IS threshold has been used since further analysis is based on ISCCP 
data. Table 1 lists these quantities for each of the 45 scenes. Some features of Fig. 2(b) are 
evident in Fig. 7(a). There is some correlation between mean LWP and cloud fraction for small 
values of the cloud fraction. There is then little discernible correlation, although overcast scenes 
tend to have the highest mean LWPs. 
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Fig. 7 (b) shows a pattern similar to Fig. 3(d) in that there appears to be little correlation 
between mean LWP and cloud fraction for partly cloudy scenes. This suggests that when scales 
of variability less than O (10 km) are averaged out, the mean liquid water path in partially cloudy 
stratus cloud fields is essentially independent of the cloud fraction. Since the total volume of 
liquid in the region is the product of mean LWP (in appropriate units) and the area covered by 
clouds, the total volume and also the mass loading per unit area of the entire region is a linear 
function of the cloud fraction. The implications of this are discussed later following an analysis of 
ISCCP data. 

4. ISCCP Analysis 

While LANDSAT data are well suited for testing the sensitivity to spatial resolution, the 
data are about 30 times as expensive to acquire and process as GOES or AVHRR data. This 
expense limited the scope of the present study to an examination of 45 oceanic cloud fields, each 
58.4-km square. More extensive analysis must rely on coarse resolution global data such as that 
cataloged in the ISCCP Cl data set (Rossow and Schiffer, 1991). 

We have restricted ourselves to one region and type of cloud in this study. In order to 
correspond most closely with the LANDSAT study, we have chosen two areas off the west coast 
of N. America during the month of July 1987. This time period coincides with the FIRE 
campaign (Cox et al., 1987) conducted in the same region. Figure 8 shows the location of the 
two 2.5° lat. x 2.5° long, grid boxes chosen for the analysis of ISCCP Cl data. 

The contours show the monthly mean planetary albedo from the Earth Radiation Budget 
Experiment scanner data (Harrison et al., 1990). The choice of areas was pardy determined by 
the availability of the most complete ISCCP Cl data on cloud properties. Other grid boxes 
tended to have missing data. Box A is at the edge of the stratus region and tended to be partly 
cloudy throughout the month while Box B was pardy cloudy or completely overcast during this 
period (Kloesel et al., 1988). 

Unlike the LANDSAT analysis, the data used here is for the same region but at three- 
hourly intervals during daylight hours. As in the previous analysis, cloud fraction at the 8 km 
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nominal resolution is one of the variables considered. An estimate of the mean radiative 
properties of the cloudy pixels is the parameter TAU cataloged in the Cl data. Details of the 
procedure used to obtain TAU are in Rossow et al. (1991). It will suffice here to mention that the 
processing is quite similar to the LANDSAT reduction described earlier, except that bi-directional 
reflectances are used in the table look up instead of the nadir reflectance given by Fig. 5. Also, 
the mean TAU for the scene is not a linear average of retrieved TAUs for each pixel but an energy 
weighted average. This is accomplished by first converting TAU to a spherical albedo before 
pixel averaging is performed. 

Figs. 9 and 10 show the scatter plot of the mean spherical albedo of the cloudy pixels in 
the scene as a function of the cloud fraction for boxes A and B respectively. There are over one 
hundred realizations presented in the plots. During the whole month, the entire scene was never 
completely clear for either of the two grid boxes. 

During each day, points represent scenes separated by only three hours temporally. In 
order to examine the degree to which cloud cover and optical properties change over that short 
span, we have identified each day of the month with a character starting with A for July 1 and 
continuing with B for July 2 and so on. The character a represents data for July 27, (3 is for July 
28, etc. With the exception of overcast days in box B, it appears that the three-hour period is 
sufficient for the cloud scene to change appreciably. 

The similarities between Figs. 9 and 10 and Fig. 3(d) are striking. Although the mean 
radiative property considered is somewhat different, both the nadir reflectance and spherical 
albedo act as proxies for the optical depth of the cloudy pixels. Unfortunately, at present, a linear 
weighted average liquid water path is not available in the ISCCP Cl data so a plot such as shown 
in Fig. 7(b) can not be presented. 

However, with the information available, the scene-averaged properties of boxes A and B 
tend to confirm the conclusions drawn from the LANDSAT analysis. At the resolution 
considered here, there is no correlation between mean cloud radiative properties and cloud 
fraction for partly cloudy scenes that are about 200 km on a side, at least for cloud fraction less 
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than about 0.8. The data for box B does show a tendency for cloudier scenes to be optically 
thicker on average when the cloud fraction exceeds 0.8. The highest mean albedos are for the 
overcast scenes. 

5. Discussion 

The study described above, though limited, suggests a direction in which cloud 
parameterizations could embark on. For low level single layer clouds, in particular maritime 
stratus, the mean liquid water path, hence optical depth, of the cloudy areas in a grid box is 
essentially invariant with cloud fraction. This appears to hold quite well for cloud fractions 
between 0.2 and 0.8. The total volume of liquid in the box is therefore a linear function of cloud 
fraction. This sort of relationship has been observed in a two-dimensional cloud ensemble model 
for the lowest model layer (Xu and Randall, 1992). Of greater relevance here, is the recent 
application of such a relationship to the modeling of surface hydrology in a GCM. 

There are some parallels between the problems involved with the parameterization of 
cloud radiative properties and surface hydrology in GCMs. Diagnostic or prognostic variables in 
GCMs are defined for the entire grid box which is typically a hundred or so kilometers on a side. 
Thus the precipitation rate or liquid water amount computed at any time step is an average value 
over the entire grid box, i.e. the expectation, E(R), of the variable R. The fate of precipitation 
falling on the surface of the model grid box, however, depends on the local precipitation rate. 
Interception, infiltration and run-off are all non-linear functions of this local rate (Entekhabi and 
Eagleson, 1989). Likewise, the radiative properties of the grid box, such as albedo and emittance 
depends crucially on the distribution of the liquid water path, not simply on the area averaged 
value (Harshvardhan and Randall, 1985). A first order approximation to the subgrid scale 
variability can be obtained by estimating the cloud fraction (or wetted area for rainfall) and mean 
cloud properties (or mean precipitation rate) of the cloudy (precipitating) area. Our study shows 
that this may be accomplished in the following manner. 

Kedem et al. (1990) and others have shown that there is a linear relationship between 
rainfall volume and fractional area of rain in convective systems. This implies that the mean 
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rainfall rate where it is raining is unique to the rainfall climatology of the location. This has been 
explained on purely statistical grounds but seems plausible. Recently, Eltahir and Bras (1993) 
have used this idea to estimate the fractional coverage of rainfall within a GCM grid box by using 
station data for the local mean rain rate and GCM output of grid-mean rain rate. 

Our stratus study has shown that there is perhaps a local mean liquid water path 
associated with partly cloudy areas of O (100 km). As with convective rainfall, it appears that at 
any instant the mean LWP of the cloudy area is the average of a population that does not change 
with change in cloud fraction, at least over the range 0.2 - 0.8. In effect, for an area the size of a 
GCM grid, there is a distribution of LWP which does not change appreciably. A simple model of 
such a cloud field is that proposed by Coakley (1991): a stratus deck is composed of cells of 
O (1 km) which tend to be optically thicker in the center and somewhat thinner at the edges. 
What appears to be a cloudy region is actually an array of such cells, each with somewhat 
different LWP distributions. For the region as a whole, the distribution could be of the form given 
by Eq. (2). Over the parameter space that this conjecture is valid, one could then estimate the 
cloud fraction based on some knowledge of the mean LWP within the cloudy portion of the box. 

There is, of course, a serious limitation to the application of this idea. Whereas local mean 
rain rate is known from station data, there is no corresponding history of cloud LWP 
observations, at least not on the global scale necessary. One possibility is to compute directly the 
mean LWP for cloudy areas based on empiricism or a cloud physics process model. After all, 
GCMs that use relationships between cloud temperature and LWP (Cess et al., 1990) are relying 
on an empirical relationship. 

The computation of a cloud fraction independently of the GCM liquid water computation 
will introduce a degree of flexibility not yet incorporated in current models. For instance, it will 
be possible to differentiate between the grid mean radiative properties of open-cell and closed-cell 
convection which are known to occur under different climatological conditions (Agee, 1987). A 
climate change experiment could then include a cloud fraction feedback which may not be trivial 
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since, for thick clouds, grid mean radiative properties are much more sensitive to cloud fraction 
than to optical depth. 
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FIGURE CAPTIONS 


Figure 1. A schematic diagram depicting the role played by threshold and resolution in 
determining the mean cloud properties of a scene. Insolation is S 0 at a particular 
zenith angle; the nadir reflectance of each pixel is I n . Increasing density of the shading 
depicts increasing reflectance. Unshaded pixels are detected to be clear. 

Figure 2. Scatter plot of the mean nadir reflectance of all the cloudy pixels in the scene vs. cloud 
fraction for LANDS AT images at a nominal resolution of 1/16 km using (a) the 
LANDSAT (LS) threshold and (b) the ISCCP (IS) threshold. 

Figure 3. As in Figure 2 but using the ISCCP (IS) threshold only for resolutions of (a) 1/16 km, 
(b) 1/2 km, (c) 2 km and (d) 8 km. 

Figure 4. Scatter plot of (a) IS - 8 vs. LS - 1/16 cloud fraction and (b) the IS - 8 vs. LS - 1/16 
mean nadir reflectance, showing the effects of using different thresholds and 
resolutions in determining the mean cloud properties of the scene. 

Figure 5. Radiative transfer results of the nadir reflectance as a function of the optical depth at 
three different solar zenith angles for plane parallel homogeneous clouds composed of 
droplets of effective radius, r e = 10 fim. 

Figure 6. Scatter plot of the inferred mean optical depth using the ISCCP threshold at 8 km 
resolution vs. the same quantity inferred using the LANDSAT threshold at 1/16 km 
resolution. 

Figure 7. Scatter plot of the mean liquid water path (mm) as a function of the cloud fraction 
derived using the ISCCP threshold (a) at 1/16 km resolution and (b) at 8 km 
resolution. 

Figure 8. The region selected for the ISCCP study. The two grid boxes chosen for the analysis 
are identified. Contours show the mean planetary albedo for July 1987 as determined 
by the Earth Radiation Budget Experiment 

Figure 9. Scatter plot of the inferred mean spherical albedo vs. the cloud fraction for daytime 
images at Box A during the month of July 1987. Characters denote the day of the 
month. See text for an explanation. 

Figure 10. As in Figure 9 but for Box B. 
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